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Abstract.  The  effective  boundary  condition  (EBC)  method  is  extended  to  nano-scale  planar  meso- 
scopic systems.  The  EBCs  appear  as  a result  of  the  2D-homogenization  procedure  andhave  the  form 
of  two-side  anisotropic  impedance  boundary  conditions  stated  on  the  structure  surface.  The  EBC 
method  supplemented  with  well-developed  mathematical  techniques  of  classical  electrodynamics 
creates  unified  basis  for  solution  of  boundary- value  problems  in  electrodynamics  of  quantum  dots. 


1.  Introduction 

The  key  peculiarities  of  quantum  dot  (QD)  heterostructures  [ ] are  related  to  spatial  con- 
finement of  the  charge  carrier  motion  and  intrinsic  spatial  in  homogeneity.  Since  the  inho- 
mogeneity scale  is  much  less  than  the  optical  wavelength,  QDs  can  be  treated  as  electrically 
small  objects  and  electromagnetic  response  of  such  structures  can  be  evaluated  by  means 
of  effective  medium  theory.  Effective  medium  approach  as  applied  to  3D  arrays  of  QDs  has 
been  developed  in  Refs.  [ , ].  In  this  paper  we  present  a general  method  which  allows  us 
to  evaluate  electromagnetic  response  of  planar  arrays  of  QDs  and  to  establish  correlation 
between  properties  of  such  systems  and  homogeneous  2D  structures  like  quantum  wells 
(QWs).  This  method,  the  effective  boundary  condition  (EBC)  method,  has  been  originally 
developed  for  microwaves  and  antenna  theory  (see,  e.g.,  [ ])  and  is  modification  of  the  ef- 
fective medium  theory  as  applied  to  2D-confined  structures.  Its  basic  idea  is  that  a smooth 
homogeneous  surface  is  considered  instead  of  the  initial  structure,  and  appropriate  EBCs 
for  the  electromagnetic  field  are  stated  for  this  surface.  These  conditions  are  chosen  in 
such  a way  that  the  spatial  structure  of  the  field  due  to  an  effective  current,  and  the  field 
of  the  real  current  in  the  initial  structure  turn  out  to  be  identical  at  some  distance  away 
from  the  surface.  Material  characteristics  of  the  structure  and  its  geometrical  parameters 
are  incl  uded  in  coefficients  of  the  EBCs. 

2.  EBCs  for  planar  nanostructures 

Under  the  derivation  of  EBCs,  the  kernel  problem  is  the  diffraction  of  electromagnetic 
field  by  an  infinite  planar  quadratic  lattice  constituted  by  identical  QDs  imbedded  in  a 
host  medium  those  permittivity  £/,  is  assumed  to  be  real-valued  and  dispersionless.  In  the 
strong  confinement  regime,  when  the  exciton  Bohr  radius  exceeds  significantly  the  QD 
linear  extention  (aB  2>  R),  the  Lorentz  dispersion  law  s(co)  — Sh  + go/io  — coq  + i/x) 
can  be  used  as  model  of  dispersion  in  a single  QD  in  the  vicinity  of  the  exciton  resonance 
[ ];  here  coq  is  the  frequency  of  the  transition,  r is  the  effective  exciton  dephasing  time  in 
the  QD,  go  = — 4t Tp.2W/HV  where  p.  and  V are  the  QD  dipole  moment  and  volume,  W is 
the  level  population  difference  (W  < 0 in  an  inverted  medium). 
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In  the  dipole  approximation,  electromagnetic  field  scattered  from  an  isolated  QD  can 
be  expressed  in  terms  of  Hertz  potentials  by 

0° 

E(r)  = Eo(r)  + £ [VV  • nfm(r)  + *?n?m(r)]  . (1) 

l,m= — oo 

Corresponding  equation  forthe  magnetic  component  is  also  presented.  Here  k = co/c.  k\  — 
k^/eh,  Eo(r)  stands  for  the  incident  field;  an  exp(— icot)  time-dependence  is  supposed. 
Further  the  incident  field  is  assumed  to  be  eA-polarized  plane  wave  propagating  at  angle  0 
with  respect  to  the  z axis.  Hertz  potentials  are  given  by 


j-re  /„■,  r exp(/*ri |r  — R/m|) 

n/m(r)  — [ eA ax x &x (R/m ) -\-  ezazz£x(Rim)j  D . . (2) 

I r — K/„,  | 

where  R/m  = {Id,  md,  0}  is  the  radius-vector  of  a lattice  site,  d is  the  lattice  period,  an 
are  the  components  of  the  QD  polarizability  tensor  a,  and  £(r)  is  the  electric  field  inside 
QD.  This  field  is  related  to  the  Hertz  potentials  by  the  equation  analogous  to  Eqs.  (1)  in 
the  limit  r — ► 0 with  the  term  l = m — 0 excluded.  For  QDs  with  planar  configuration  in 
the  vOy  plane,  e.g.,  discs,  islands,  flattened  pyramids,  etc.,  the  QD  polarizability  in  the  z 
direction  azz  can  be  neglected. 

The  next  step  is  the  2D-averaging  of  the  electromagnetic  field  in  the  z = 0 plane.  Such 
a procedure  implies  replacement  of  the  discrete  2D  elementary  scatterer  by  a homogeneous 
d x d surface  element;  mathematically  this  procedure  reduces  summation  in  equations  for 
inner  and  outer  fields  to  integration.  Then,  one  should  find  discontinuities  of  the  mean  field 
tangential  components  atz  = 0 taking  into  account  that  the  mean  field  has  exp(ik\x  sin  6) 
dependence  on  the  x co-ordinate.  Then,  after  some  manipulations,  we  come  to  the  covariant 
notation  of  the  EBCs: 


n x n x (H1  - H11)  = x (E1  +En), 

n x (E1  — E11)  = -fnx  V[n  • (E1  + E11)] 


(3) 


where 


er  = 


Here,  I y is  the  2 x 2 unit  tensor  and  «||  is  the  surface  polarizability  tensor  with  components 
atj  ( i , j — x,  y);  Sx  « — 8/>/2 d , Sz  as  lit yfn /d  . 

The  equations  (3)  and  (4)  constitute  the  complete  system  of  EBCs  for  electromagnetic 
field  in  planar  QD  structures.  The  technique  of  macroscopic  averaging  used  under  their 
derivation  is  similar  to  one  which  introduces  the  constitutive  parameters  for  bulk  media, 
but  differs  in  that  the  averaging  occurs  in  boundary  conditions  rather  than  in  field  equa- 
tions. Thus,  in  electrodynamics  of  low-dimensional  struct  ures  the  EBCs  play  the  same  role 
as  constitutive  relations  in  electrodynamics  of  bulk  media.  The  EBCs  keep  validity  for 
arbitrary  configuration  of  elementary  cell  and  for  planar  layers  with  random  distribution 
of  QDs.  The  difference  will  manifest  itself  in  the  modified  coefficients  <5,-.  Since  EBCs 
(3)  are  local  and  the  coefficients  <X||,  § do  not  depend  on  the  angle  of  incidence  0,  they 
are  applicable  at  arbitrary  excitation  in  spite  of  that  they  have  been  obtained  originally  for 
plane  waves.  Moreover,  the  EBCs  (3)  are  analogous  to  corresponding  EBCs  for  QWs  (see, 
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e.g.,  [ ] ) if  spatial  dispersion  in  the  latter  can  be  neglected.  Thus,  a planar  layer  comprising 
a 2D  array  of  QDs  can  be  treated  as  an  effective  Q W.  As  a result,  well-developed  formalism 
of  investigation  of  QWs  can  be  extended  to  QD  arrays  by  introducing  of  effective  integral 
parameters  of  the  array  defined  by  Eqs.  (4). 

3.  The  role  of  nonloeality 

In  the  weak  confinement  regime,  the  QD  electromagnetic  response  becomes  nonlocal  and 
the  medi  um  polarization  takes  the  form  of  the  integral  operator  [ ]: 

P(r)  = A<t>(r)  J dKrQEtr')  J V , (5) 

v 

where  A ~ 1 /(co—coo+i/t ).  <J>(r)  is  the  envelope  function  of  the  exciton  ground  state.  Eq. 
(5)  defines  very  special  type  of  nonlocality:  the  integral  operator  has  degenerated  kernel. 
This  property  makes  possible  analytical  consideration  of  the  nonlocality  problem:  integral 
differential  equations  describing  electromagnetic  field  in  QD  turn  out  to  be  equivalent 
to  the  integral  Fredholm  equations  with  degenerated  kernels.  As  a result,  we  obtain  the 
golarizability  tensor  of  an  isolated  QD  as  S = AN2  (l  — 4tt  A Y)“’  , where  the  3D-tensor 
Y is  given  by 

1 C C d2  I 

Tap  = — / <t>(r)<t>(r')— — — -J3rdV  xa,xp=x,y,z.  (6) 

4 n J J dxadxp  |r  - r'| 

v v 

Thus,  we  have  shown  that  the  nonlocality  changes  values  of  the  polarizability  tensor  compo- 
nents but  does  not  change  the  general  representation  of  the  scattering  operators  as  compared 
to  the  strong  confinement  regime.  This  allows  us  to  conclude  that  the  above  introduced 
EBCs  remain  valid  in  the  weak  confinement  regime  as  well.  Note  that  the  above  result 
admits  extension  of  the  Maxwell  Garnett  approach  to  3D  composites  constituted  by  QDs 
in  weak  confinement  regime. 


4.  Radiative  decay  rate  in  planar  array  of  QDs 

Let  us  apply  the  EBC  method  for  the  investigation  of  exciton  radiative  time  in  2D-array 
of  QDs  which  are  assumed  to  be  spherical  inclusions  of  the  radius  R.  The  problem  of  the 
radiative  linewidth  evaluation  in  QWs  is  solved  by  finding  of  the  frequency  poles  of  the 
reflection  coefficients  [ ].  Real  parts  of  these  poles  determine  resonant  frequencies  while 
imaginary  parts  give  the  homogeneous  linewidths.  It  can  easily  be  shown  that  EBCs  (3) 
with  <X||  = ickr]l\\/4it,  § = —r]/2s(co)  and  tj  — Lqw[s(co)  — €/,]  describe  a QW  with 
Lq w as  its  thickness.  Thus,  the  reflection  coefficient  for  planar  array  of  spherical  QDs  are 
given  by  a corresponding  equation  for  QW  after  substitutions  Lqw  —*■  2 R,  g®W  —*■  gQff  « 
(2tt  R2 /3d2)go  and  coo  —*■  &>o  — coo  — go/3sh-  Then,  using  results  of  Ref.  [ ] and  above 
defined  substitutions,  we  obtain  the  expressions  r/;  = To/ cos  9 , rv  = To  cos  9 , for  T E 
and  T M polaritons,  correspondingly;  here 


r0« 


2tx  R2  coo 
Sy/eJ,  d 2 c 


Rgo , 


(7) 


Comparison  of  T o with  the  radiative  decay  rate  of  a single  spherical  QD  y gives  us  Tq/y  — 
B « 3n/(k\d)2.  In  dense  arrays  of  QDs  B 1 signaling  significant  enhancement  of 
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the  decay  rate  Analogous  superradiance  factor  with  the  Bohr  radius  as  instead  of  d was 
introduced  also  for  QWs  [ ].  One  can  interpret  the  coefficient  B in  close  analogy  with  Ref. 
[ ]:  it  results  from  the  coherent  excitation  of  QDs  located  at  distance  d from  each  other. 

5.  Conclusion 

EBCs  given  by  Eqs.  (3)  state  mathematical  equivalence  of  optical  properties  of  a 2D  pe- 
riodical layer  of  QDs  and  an  isolated  QW.  It  should  be  stressed  that  the  mechanisms  of 
transport  processes  and  oscillator  strengths  in  each  case  are  essentially  different.  Never- 
theless, the  equivalence  makes  it  possible  to  extend  to  QD-based  planar  structures  with 
more  complicated  configurations  (finite-sized  QD  layer,  QD  layer  in  microcavity,  several 
QD  layers,  etc.)  the  well-developed  mathematical  formalism  of  investigation  of  QWs. 
Namely  this  equivalence  provides  promising  potentiality  of  the  derived  EBCs  for  particu- 
lar electrodynamical  problems  in  QD-based  structures.  In  particular,  threshold  current  for 
QD-based  lasers  can  be  evaluated  by  analogy  with  solution  of  corresponding  problem  for 
the  QW  lasers;  the  EBC  method  allows  us  to  analyze  electromagnetic  response  of  a QD 
layer  (or  a multilayer  structure)  placed  in  microcavity:  this  is  very  important  for  the  design 
of  QD-based  semiconductor  lasers  [ ]. 
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